1 


"C 


CO 

LI 

o 

£ 


< 

O 

t 


/vTty, 

& 




f£B 






** — 


m A'frt 187 t 1 

£ fr^y-No v 2 4 

RM No. L8T11 


NACA 



(jjt 

\r-KJi KO^O >- 

R„U.(.v-U > 5 

1^ lW 


K: ! » 

‘CTED 



RESEARCH MEMORANDUM 


INITIAL EXPERIMENTS ON FLUTTER OF UNSWEPT 
CANTILEVER WINGS AT MACH NUMBER 1.3 

By 


ler 


W. J. Tuovila, John E. Baker, and Arthuj^Jr^^ 

Langley Aeronautical Laborator^£ ^ ^ 




to 


Srcr 


V7F — Authority 

'a* oij Sfco cJ 


/3,/fS-H 


Date 

% 4 1354 




Langley Field, Va. 


CLASSIFIED DOCUMENT 


This document contains classified Information 
affecting the National Defense of the United 
States within the meaning of the Espionage Act, 
USC 50:31 and 32. Its transmission or the 
revelation of its contents in any manner to an 
unauthorized person is prohibited by law. 

Information so classified may be imparted 
only to persons in the military and naval 
services of the United States, appropriate 
civilian officers and employees of the Federal 
government who have a legitimate Interest 
therein, and to United States citizens of known 
loyalty and discretion who of necessity must be 
informed thereof. 


v 

* v 


y 




*/ / 


TV 






% 




NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON jan 1 1 1949 


January 6, 1949 




l 


a 


-umfissF© 


UNCLASSIFIED 



NACA RM No. L8jll 


UN&AS8EE 

IT 



Clasnilca 


NATIONAL ADVISORY COMMITTEE FOR A 


RESEARCH MEMORANDUM 



RLCTED 

jd-i, |-j, 


J 


‘UJ 




INITIAL EXPERIMENTS ON FLUTTER OF UNSWEPT 
CANTILEVER WINGS AT MACH NUMBER 1.3 
By W. J. Tuovila, John E. Baker, and Arthur A. Regier 


SUMMARY 


A supersonic tunnel designed to operate at Mach number 1.3 was 
used for a preliminary experimental flutter investigation of widely 
different unswept cantilever wings. Data for 12 wings with mass- 
density parameters 1 /k ranging from 52 to 268 , center— of-gravity 
positions ranging from 46 to 63 percent chord from the leading edge, 
and elastic— axis positions ranging from 3 ^ to 52 percent chord from 
the leading edge are considered. 

A comparison is made of the test results with calculations of 
"bending— torsion flutter obtained by the theory of flutter in supersonic 
two-dimensional flow and it is concluded that the test data are in 
reasonable agreement with the calculated results. In general, the 
theoretical values are conservative. As shown by the theory, the flutter 
results are quite sensitive to the location of the center of gravity. 
Thick and thin, blunt and sharp airfoil-section shapes were used, but 
no very pronounced effect of the section shape on flutter characteristics 
was found. The data suggest that for cantilever wings the bending degree 
of freedom may suppress the one— degree— of— freedom torsional flutter and 
that coupled bending— torsion flutter effects occur. The experiments 
include a study of the effect of the addition of tip moments of inertia. 
With the center of gravity of the tip weights coincident with the center 
of gravity of the wing section no detrimental effect on the flutter 
speed was found. 
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P 

M 

1/k 

*a 

r a 2 

V 

f f 

f h 


mass density of air in test section 
mass of wing per unit span 
mass— density parameter 

mass moment of inertia of wing about elastic axis per unit 
span 


(v^ 2 ) 


radius— of-gyration parameter 

flutter velocity, feet per second 
flutter frequency, cycles per second 

first bending frequency, cycles per second 
first torsion frequency, cycles per second 


o>h = 2*f h 


m = 2nf rf 
a a 


Sh 

6 a 


first bending damping coefficient 
first torsion damping coefficient 


INTRODUCTION 


The background and theory for the flutter of an airfoil in a two- 
dimensional flow at supersonic speeds is given in reference 1 . The 
present investigation is a preliminary survey to determine the possibility 
of using the theory of reference 1 for flutter at supersonic speeds to 
predict the coupled bending— torsion flutter of widely different unswept 
cantilever wings at a low supersonic Mach number. This preliminary 
investigation is not intended as a critical test of the theory since the 
analysis does not consider the effect of mode shape, aspect ratio, 
section shape, tip Mach cone, or viscous effects. 

A single— degree-of— freedom torsional instability which may occur in 
the Mach number range 1.0 to 1.58 is discussed in reference 1 . In order 
to also investigate the possible occurrence of such single— degree flutter 
on cantilever wings, the test apparatus was designed to operate at a 
Mach number of 1.3. 


CONFIDENTIAL • 



NACA RM No. l8jll 


3 



A brief discussion is given of the effects of concentrated masses 
placed at the wing tip, the center of gravity of the masses coinciding 
with the center of gravity of the wing, and the effects of sharp and 
blunt leading edges on the wings. 


APPARATUS AND TEST METHODS 

The tests were made at a Mach number of 1.3 in an "intermittent" 
two-dimensional supersonic tunnel having a 9*24— inch by 18. 23— inch 
test section (figs. 1 and 2). The apparatus operates from atmospheric 
pressure to a vacuum. A quick-operating valve allows a steady— flow 
condition to be reached in approximately 0.15 second and this condition 
of steady flow can be maintained for 5 ‘to 7 seconds. To prevent con- 
densation in the test section, a room was constructed at the tunnel 
entrance in which the air could be heated. Variation of the air 
temperature from 170° F to 220° F results in a velocity range at the 
test section from 950 miles per hour to 990 miles per hour at a Mach 
number of 1.3. The test-section Mach number determined by optical means 
varied from 1.29 to I.3I. The test— section Mach number determined by a 
pressure survey showed a variation from 1.31 to 1.33 (tig* 3)* 

The models were mounted cantilever fashion from the side of the 
tunnel (fig. 4 -). In order to avoid possible model failure during the 
tunnel transient conditions, and since supersonic flutter characteristics 
were being determined, it was desirable to withhold the model from the 
flow during the periods of acceleration and deceleration. A pneumatic- 
cylinder arrangement was installed for this purpose (fig. 2). By using 
this device the model could be held outside the tunnel wall until stable 
flow conditions were reached. Then, by means of electrically operated 
valves, controlled by an observer, the model could be allowed to enter 
the tunnel slowly and to withdraw quickly if necessary. 

The flutter models having rectangular plan forms were constructed 
of laminated Sitka spruce and also of duralumin. The wing dimensions 
and parameters are given in table I. The chords ranged from 3*03 to 
4.22 inches and the lengths ranged from 6 to 9*12 inches. Both thick 
and thin sections with blunt and sharp leading edges were used. The 
airfoil sections used were 3-, 5-, and 8-percent-thick circular arcs, 

3— percent— thick double wedge, NACA l6— 010, and NACA 65— 007. The mass- 
density parameter 1 /k ranged from 52 to 268, the center of gravity 
ranged from 46 to 63 percent chord, and the elastic axis ranged from 34 
to 52 percent chord. 

Before each model was installed in the tunnel, its weight, moment of 
inertia, and section center-of-gravity position were determined. After 
installation in the tunnel the elastic axis was located by observing the 
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axis of zero twist optically. The first tending frequency and the 
damping were obtained from a free— vibration record of the model. Since 
the wings were uniform, the stiffness parameter GJ could be computed 
from a torsional— vibration record obtained with a mass of large known 
moment of inertia attached to the wing tip. The uncoupled first torsion 
frequency could then be calculated. The torsional damping was determined 
from the torsional— vibration— decay curve. 

The models were designed not to flutter on the first run. The 

models were progressively modified by adding lead to shift their centers 

of gravity and by slotting to shift their elastic axes and change their 

bending and torsion frequencies, until flutter occurred. If a model was 

saved, its parameters were changed to yield another flutter point. The 

values of the radius— of-gyration parameter r ^ and center of gravity 

cc 

were determined from the unmodified wing and the added masses. 

During each test run, the following data were recorded simultaneously 
by means of a recording oscillograph: flutter frequency, position of the 

model in the tunnel, and static pressure which indicates the Mach number. 

A sample record of the flutter of model B— 5 is given in figure 5 . 


METHOD OF ANALYSIS 


The flutter data for the particular wings tested are put in coef- 
ficient form and compared with the analytic solution of the supersonic 
bending— torsion flutter problem in a two-dimensional flow given in 
reference 1. The data of this paper were obtained at a Mach number 
of 1.31 and, since aerodynamic coefficients at this Mach number are not 
included in reference 1 , a linear interpolation was made between the 
computed values of the flutter— speed coefficient at Mach numbers of 1.25 
and 1.^3, for which the aerodynamic coefficients are available. First 
bending and uncoupled first torsion frequencies and damping 
coefficients g h and g a were used in the theoretical analyses. 

It is desirable to examine the possible errors introduced into the 
results by the method of interpolation. It is known that the error may 
be very large; for example, in the single— degree— of— torsional— instability 
case for the elastic-axis position at 50 percent chord the interpolation 
was made directly between the aerodynamic coefficients at Mach 
numbers 1.25 and 1.43* This was necessary since the wing is stable at 
a Mach number of 1.43 an & the theory yields no solution for the flutter- 
speed coefficient under these conditions. 
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RESULTS AND DISCUSSION 


The significant flutter parameters and the results of the calcu- 
lations are given in table I. As indicated by test methods, a large 
number of runs were made on wings which did not flutter. Table I gives 
only the results for the wings for which flutter was observed. In 
altering the model to obtain flutter, the center of gravity was moved 
back in steps of about 2 percent of the chord. Consequently, the 
position of the center of gravity which would produce flutter is known 
to about 2 percent of the chord and lies between the position that did 
not lead to flutter and the position that produced flutter. The results 
are sensitive to center— of-gravity position and this may account for 
some of the scatter of the data. Contributing also to the scatter of 
the data are the inaccuracies in obtaining the wing parameters, the 
effect of the degree of penetration of the model into the tunnel, and 
errors in the determination of the flutter— speed coefficients arising 
from the method of interpolation. 

Some of the results listed in table I are presented in figures 6 
and 7- In figure 6 the theoretical and experimental flutter— speed 
coefficients are compared. The fact that the experimental data fall 
above the k5° line, for the most part, indicates that the theory of 
reference 1 is generally conservative as far as application to cantilever 
wings is concerned. Considering that a slight inaccuracy in the location 
of the center of gravity has a large effect and considering also unac- 
counted effects of section shape, aspect ratio, and Mach cone, the 
agreement is not unsatisfactory. The theoretical and experimental 
flutter frequencies are compared in figure 7 ; the experimental frequencies 
ranged from about 0.6l to 1.09 times the theoretical values. In all cases 
the mode at flutter appeared to consist mainly of a coupling of first 
bending and first torsion modes. 

Since the present investigation is of a preliminary nature and 
covers a wide range of parameters, no attempt was made to isolate the 
effects of separate parameters such as the mass-density parameter, 
frequency ratio, elastic axis, and center of gravity which are treated 
by the two-dimensional theory, or parameters such as aspect ratio not 
treated by the theory. 


An attempt was made to investigate some of the possible effects of 
airfoil section shape on flutter. The intermingling of the data for the 
various models (figs. 6 and 7) suggests that the section shape has no 
very pronounced effect on flutter at Mach number 1.3* However, more 
difficulty due to divergence was experienced with thick models having 
blunt leading edges. This is in accord with the higher— order method of 
calculation (order higher than in the linear method) for pressure distri- 
bution at supersonic speeds in steady flow which shows that the center 
of pressure may move ahead of the 50-percent-chord position for thick 
blunt-nosed airfoils particularly at Mach numbers near unity. It was 
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observed in the tests that the thick airfoils tended to diverge even 
though the elastic axis was ahead of the 50 - percent-chord position. 

Since practical winged vehicles pass through the subsonic speed 
range in order to reach supersonic speeds, some discussion of and 
comparison with subsonic criterions are desirable. For this purpose, 
incompressible flutter— speed coefficients were computed by the method 
of reference 2 for which first bending and uncoupled first torsion 
frequencies and damping coefficients g^ and g a were utilized. That 

flutter— speed coefficients based on subsonic theory are conservative 
with respect to supersonic speeds has been suggested in reference 3 and, 
in fact, is indicated by some of the calculations in reference 1. This 
is also indicated by the present tests, as illustrated in figure 8 in 
which the experimental flutter-speed coefficients are plotted against 
the incompressible theoretical flutter— speed coefficients. The statement 
may not be true in general; for example, the condition when the frequency 

ratio — =1 may need further investigation and, in any case, the margin 

of safety is not established. Some of the models were permitted to 
encounter the tunnel transient speeds and, for example, model "F" which 
had fluttered at Mach number 1.3 was held in the tunnel while the tunnel 
was brought up to speed. The wing fluttered and broke at a Mach number 
of about 0.7, a result which is in general agreement with the subsonic 
calculation. Flutter data obtained with bombs and rocket missiles 
(references 4 to 6) and other experience indicate that if flutter failures 
occur, they occur, in general, at near sonic speeds. For the practical 
purpose of making preliminary estimates of a wing flutter speed such 
formulas as, for example, the approximate flutter formula in reference 2, 
or similar criterions, thus appear useful over a wide range of speeds. 

In reference 3, Smilg suggests a torsional frequency crite- 
rion o^Cy- > 1000 feet per second as sufficient to prevent one degree of 
torsional flutter at supersonic speeds. The criterion is based on the 
assumption that the first bending frequency is very high with respect to the 
first torsion frequency. In order to look into this criterion a cantilever 
model was equipped with tip weights at both the leading and trailing edges 
which reduced the torsional frequency. The results of the tests are 
shown in table II. In all cases a slight shift of the center of gravity 
ahead of the. location at flutter stopped the flutter. The fact that 
flutter is extremely sensitive to the center— of-gravity position and 
that the values of the product u) a c w are far below 1000 indicates that 

the criterion is overly conservative when applied to cantilever wings 
with normal bending— torsion frequency ratios. The data further suggest 
that for cantilever wings the bending degree of freedom may suppress 
the one— degree— of— freedom torsional flutter and that bending— torsion 
effects occur. 
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The data further indicate that no harmful effect on the flutter 
speed ensues when the center of gravity of the tip weights and the wing 
coincide. It is observed that the frequency ratio varies from 0.55 to O.85 
and that the torsional frequency has been reduced to as low as one— third 
the bare wing value. For the largest mass moment of inertia on the wing 
tip (last case in table II) it was necessary to move the center of gravity 
farther toward the trailing edge to produce flutter. 

In figure 9 * the curves represent theoretical flutter— speed coef- 
ficients for one— degree— of— freedom flutter according to reference 1. 

For the sake of possible interest, the experimental flutter— speed coef- 
ficients corresponding, however, to the coupled bending— torsion values 
are shown superposed on the figure. 

It was hoped that some systematic aspect— ratio effects could be 
obtained from the present tests, but the results were rather contradictory. 


Some models were used which spanned the tunnel ( except for inch tip 


clearance j so that presumably two-dimensional flow over the wing could 

be expected. By retracting the tip from the boundary layer, flutter of 
full— span models of NACA 16-010 section could be stopped; however, by 
retracting the tip from the boundary layer, the flutter amplitude of 
^-percent double-^wedge models increased. The effect of the subsonic 
boundary layer at the tip of the models is not known. In one particular 
case the model fluttered at 99 percent of the theoretical frequency on 
entering the tunnel and the frequency changed to 128 percent at a smaller 
amplitude when the model spanned the tunnel. As the model was retracted, 
the flutter frequency dropped to 99 percent of the theoretical value and 
fluttered to destruction. A more systematic investigation of the aspect 
ratio, and tip and shape effects is desirable to clarify various aspects 
of the problem. 


1 . There is rather satisfactory agreement between experimental and 
calculated flutter-speed coefficients. In general, the theoretical values 
are conservative. 

2 . No very pronounced effect of airfoil— section shape on the flutter 
characteristics was found in these preliminary experiments; however, 
significant divergence effects were observed on thick blunt-nosed airfoils. 

3. The data suggest that for cantilever wings the bending degree of 
freedom may suppress the one degree of torsional flutter and that coupled 
bending-torsion effects occur. 
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4. More systematic experimental examination of the effects of the 
individual parameters is desirable. 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I 


EXPERIMENTAL AND THEORETICAL RESULTS OF FLUTTER INVESTIGATIONS 


Model designation 

Parameters ' 

A-l 

B— 1 

B— 2 

B— 3 

B— 4 

B-5 

C-l 

C— 2 

D— 1 

E— 1 

F— 1 

G— 1 

Description of model section . . . 

NACA 

65-007 

NACA 

16-010 

NACA 

16-010 

NACA 

16-010 

NACA 

16-010 

NACA 

16-010 

8 percent 
circular arc 

8 percent 
circular arc 

Modified 
circular arc 
4.74 percent thick 

3 percent 
circular arc 

3 percent 
double wedge 

Modified 
circular arc 
5 percent 
thick 

Length, in 

6.06 

9.125 

9.125 

9.125 

9.125 

7.50 

6.00 

6.00 

7.125 

9.125 

9.125 

6.00 

Chord, in 

3.03 

4.03 

4.03 

4.03 

4.03 

4.03 

3.03 

3.03 

4.22 

4.04 

4.01 

4.00 

Center of gravity, percent chord . 

1*9.1 

51.6 

54.6 

56.0 

56.7 

57-0 

. 53.0 

55.4 

46.0 

62.6 

56.5 

53-5 

Elastic axis, percent chord . . . 

41.3 

34.1 

39.6 

39.6 

44.2 

39*5 

48.0 

51.55 

37.0 

38.7 

45.2 

47.5 

k a 

0.28 

0.39 

O .38 

0.40 

0.37 

0.37 

0.230 

0.233 

0.275 

0.510 

0.29 

0.27 

1 /k 

64.9 

95.3 

108.1 

113.1 

113.3 

130 

67.1 

74.1 

53.5 

267.5 

150.8 

51.7 

GJ 

578 

3060 

2945 

3150 

4105 

2500 

839 

909.5 

1025 

2710 

2220 

1066 

V cp8 

133.4 

95.0 

100.9 

98.0 

104.2 

104.0 

157 

157 

110 

25.2 

28.5 

148.5 

cps 

278 

I 63 

156.3 

154.8 

183 

162 

363 

361 

178 

81.7 

132.8 

245 

ff, cps 

180 

134 

132.5 

134.6 

131.5 

136 

188 

184 

142 

71 

70.5 

176 

p X 10^ (test section) 

0.876 

0.945 

0.890 

0.897 

0.908 

0.886 

0.887 

0.884 

0.918 

0.888 

0.888 

0.910 

g h 

0.04 

0.03 

0.025 

0.03 

0.025 

0.035 

0.05 

O .05 

0.03 

0.004 

0.01 

0.04 

g a 

0.04 

0.04 

0.04 

0.035 

0.03 

O.O 35 

0.04 

0.04 

0.03 

0.005 

0.005 

0.04 

“h/“a • 

0.48 

0.583 

0.645 

0.633 

0.57 

0.64 

0.432 

0.435 

0.614 

O .308 

0.215 

0.606 

V/bajj. (experimental) . . . f . . . 

10.15 

9.98 

IO .31 

10.20 

10.40 

IO .03 

9.715 

9.92 

9.04 

19.13 

19.61 

7.71 

V/boif (theoretical) 

6.946 

7.770 

7.109 

7.546 

7.938 

7.683 

8.53 

7.18 

5.28 

14.85 

9.38 

4.92 

V/bo> a (experimental) * . 

6.59 

8.21 

8.74 

8.91 

7.44 

8.42 

5.02 

5.055 

7.25 

16.7 

10.35 

5.54 

V/bao a (theoretical) ....... 

6.554 

7.158 

5.664 

6.029 

6.234 

6.041 

4.913 

5.503 

5.069 

13.08 

8.23 

3.905 

(o^/co^ (experimental) ....... 

0.648 

0.822 

0.847 

0.870 

0.719 

0.840 

0.518 

0.51 

0.798 

0.868 

0.531 

O .718 

o>f/a> a (theoretical) 

0.944 

0.921 

0.797 

0.800 

0.785 

0.786 

0.577 

O .766 

O .96 

0.88 

0.868 

0.795 

(V/bo) a ) (M » 0) (theoretical) . . . 

3.427 

4.302 

4.245 

4.358 

4.121 

4.612 

3 . 184* 

3-374 

2.811 

8.263 

5.491 

2.585 
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TABLE II 


FLUTTER DATA ON A CANTILEVER MODEL WITH TIP WEIGHTS 


1 2 
Wing chord = — ft; elastic axis = 47 percent from L. E.; r a = 0.23 (for hare wing 

L 3 

weights); weight of hare wing = 0.0806 lh; length of wing = 6^ in.; tip-weight 
gravity coincides with the wing center of gravity 


with no tip 
center of 


I — ■ — ' 

Frequency 
( cps ) 

Center of gravity 
of wing and 
tip weights 
(percent chord) 

c w CD a 

Tip weights 
(lh) 

Moment of 
inertia of 
tip weights 
about c. g. 

(in. -lb-sec ) 

Torsion 

Bending 

Flutter 

w 

• 

T. E. 

244 

134 

153 

50 

510 

0 

0 

0 

133 

102 

105 

50 

• 

278 

.00949 

.00949 

.000162 

103 

81 

86 

. 49 

216 

.01766 

.01754 

.000342 

80 

68 

74 

53 

167 

.02747 

.03044 

.000686 
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Figure 1 . - Diagram of the supersonic flutter apparatus . 
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Figure 2 . - 


General view of the Langley supersonic flutter apparatus . 
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Figure 3*" Mach number survey of the test section. 
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Figure 4.- Flutter model installed in the test section. 
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Figure 5.- Sample oscillograph record of the flutter of model B -5. £ 
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Figure 6 . - Comparison of experimental values of V to theoretical 

values of V /ha) a at M = 1.31* 
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Figure 7 . - Comparison of experimental 

values of 


values of 
at M = 1.31. 
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Figure 8.- Comparison of experimental values of V /txn a at M = 1*3 
with theoretical incompressible values of V /b<D a at M = 0. 


(JNoLrtC>6lri£i?* 



Figure 9.- Comparison of experimental values of V/h for coupled wing 
flutter with theoretical values for one -degree- of -freedom torsional 
flutter. M = 1.3, 80 , = 0.0k. 
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